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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 

UNITS  OF  MEASUREMENTS 

U.S.  Customary  Units  of  Measurement  Used  in  this  Report  can  be 
converted  to  metric  (SI)  units  as  follows: 


Multiply 


By 


To  Obtain 


i  nches 
feet 

cubic  feet  per  second 
pounds  (mass) 

pounds  (mass)  per  cubic  foot 


25.400 

0.3048 

0.0283 

0.4535924 

0.15709 


mi  1 1 imeters 
metres 

cubic  metres  per  second 
kilograms 

kilonewtons  per  cubic  metre 
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STUDY  OF  EMBANKMENT  PERFORMANCE  DURING 
OVERTOPPING  AND  THROUGHFLOW 


Phase  I:  Rigid  Model  Studies 
Phase  II:  Erodable  Model  Studies 

PART  I:  INTRODUCTION 

1.  Overtopping  of  embankments  and  earth  dams  has  been  a  major 
cause  of  failure  and  has  also  resulted  in  millions  of  dollars  of  damage 
in  recent  years.  Data  on  actual  overtopping  events  is  scarce,  and  full 
scale  tests  to  study  the  phenomenon  are  very  expensive.  Because 
gravity  induced  factors  govern  the  failure  mechanism,  centrifgual 
modeling  was  proposed  to  study  overtopping,  and  the  feasibility  of  this 
method  was  shown  in  the  initial  phase  of  this  study  (Ko  et  al.,  1984). 

2.  In  the  final  stages  of  the  feasibility  study  questions  were 
raised  regarding  the  appropriateness  of  centrifuge  modeling  of  hydrau¬ 
lic  phenomena  such  as  open  channel  flow.  Therefore,  it  was  decided  to 
perform  additional  centrifuge  tests  as  a  modeling  of  models  study  using 
rigid  non-erodable  embankment  models  to  allow  for  the  measurement  of 
hydraulic  parameters,  notably  flow  quantity  and  velocity,  in  the 
centrifuge  and  to  verify  the  scaling  relations  for  these  parameters. 

3.  The  validity  of  any  modeling  experiments  is  best  demonstrated 
by  comparison  with  full  scale  prototype  performance,  and  at  the  time 
the  additional  centrifuge  tests  were  being  planned,  Simons,  Li,  and 
Associates  (SLA)  of  Fort  Collins,  Colorado  were  constructing  a  series 
of  full  scale  tests  of  highway  embankment  sections  in  a  study  of 
overtopping  erosion  for  the  Federal  Highway  Administration  (FHWA). 

Thus  it  was  decided  to  model  one  of  the  SLA  test  sections  in  the 
centrifuge  tests  in  a  two-phase  test  program.  Phase  I  consisted  of  the 
non-erodable  embankment  tests,  while  Phase  II  utilized  the  same 
embankment  geometry  in  erodable  embankment  models,  constructed  of  the 
sandy  clay  soil  used  to  construct  the  SLA  embankments,  and  provided 
data  on  erosion  rates  for  comparison  to  prototype  data. 
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PART  II:  EQUIPMENT  AND  INSTRUMENTATION 


The  Centrifuge 


4.  The  10  g-ton  geotechnical  centrifuge  which  has  been  in 
operation  at  the  University  of  Colorado  at  Boulder  since  1978  is  shown 
schematically  in  Figure  1 1- 1 .  This  centrifuge  is  a  Genisco  Model 
1230-5  G  accelerator  with  a  rebuilt  symmetrical  rotor  to  carry  swinging 
baskets  at  the  ends.  The  radius  to  hinges  where  the  baskets  are 
attached  is  42  in.  and  that  to  the  bottom  of  the  baskets  in  the 
swung-up  positions  is  53.5  in.  One  basket  carries  the  test  payload, 
while  the  other  carries  the  counterweights.  The  maximum  payload  size 

i  s  18  i  n.  x  18  i  n.  x  1?  i n.  The  rated  capacity  of  the  centrifuge  is  10 
g-tons,  typically  100  g  on  a  200  lb  payload. 

5.  The  centrifuge  is  equipped  with  56  electrical  slip  rings  and 
two  hydraulic  rotary  joints  which  have  0.250  in.  ID  tubing  inside  the 
centrifuge.  As  these  rotary  joints  were  judged  inadequate  in  supplying 
the  flow  rates  required  for  the  overtopping  experiments  in  this 
project,  an  additional  large  capacity  toroid  water  conveyance  system 
was  added  to  the  centrifuge  as  part  of  the  feasibility  study.  This 
system  is  described  in  the  original  report  (Ko  et  al.,  1984). 

6.  The  system  was  modified  for  the  tests  described  herein  to 
allow  for  measurement  of  flow  quantity  and  improved  control  of  flow  by 
using  two  shut  off  valves  as  shown  in  Figure  1 1-2.  Both  valves  1  and  2 
were  opened,  with  the  end  cap  closed,  to  direct  flow  to  the  centri¬ 
fuge.  Valve  1  was  used  for  flow  quantity  control  and  thus  overtopping 
depth,  while  valve  2  was  fully  opened.  To  measure  flow  quantities 
after  an  overtopping  test,  valve  2  was  closed  while  valve  1  was 
untouched,  and  the  end  cap  removed.  The  flow  quantity,  Q  in  cfs,  was 
measured  by  collecting  the  flow  in  a  30  gallon  tank  over  a  time  period, 
usually  approximately  20  seconds,  and  weighing  the  collected  water. 


-5- 


Sample  Container 


7.  The  hydraulic  modeling  and  erosion  experiments  were  conducted 
in  a  new  sample  container  constructed  specifically  for  these  tests. 

This  container  is  shown  schematically  in  Figures  1 1-3  and  1 1-4.  Due  to 
the  limited  space  in  the  centrifuge  swinging  basket,  only  a  portion  of 
the  SLA  embankment  crest  and  downstream  slope  could  be  modeled  in  most 
of  the  tests.  One  erodable  embankment  test  was  completed  using  the 
full  SLA  embankment  configuration  and  is  described  in  a  later  section. 
To  accommodate  models  of  various  widths  used  in  the  modeling  of  models 
study,  one  wall  of  the  sample  container  was  adjustable. 

8.  As  in  the  feasibility  study,  the  water  was  stored  in  a 
storage  reservoir  located  behind  the  embankment  section  and  after 
overtopping,  the  water  drained  through  drain  ports  in  the  compartment 
below  the  support  floor.  The  crest  of  each  test  embankment  was  at  the 
same  elevation.  The  water  was  transmitted  from  the  toroid  catchment 
through  a  1.25  inch  O.D.  PVC  pipe  with  a  diffuser  at  its  end,  to 
minimize  turbulence  of  water  as  it  flowed  over  the  crest.  To  allow 
observation  of  the  flow  cross  section  and  embankment  erosion,  a  1.0 
inch  thick  plexiglass  plate  was  used  as  a  side  wall  of  the  sample  con¬ 
tainer  as  shown  in  Figure  1 1-4.  A  one-half  inch  square  grid  pattern 
was  inscribed  on  this  plate  to  provide  a  reference  for  measurement  of 
flow  depth  and  erosion  depth  in  the  photographs  taken  during  the 
tests.  Various  basket  mounting  positions  in  the  centrifuge  were  exam¬ 
ined  and  the  one  adopted  for  these  tests  is  shown  in  Figures  1 1 - 3  and 

1 1-4.  Mounting  the  container  in  any  other  position  is  not  desirable 
since  the  free  water  surface  would  be  curved  due  to  a  divergent  gravity 
field,  hence  making  it  difficult  to  read  the  overtopping  depth.  As 
shown  in  Figure  1 1 -5,  placing  the  models  parallel  to  the  artificial 
gravity  field  creates  a  smaller  angle,  y,  of  divergent  accelerations 
compared  to  placing  the  models  perpendicular  to  the  artificial  gravity 
fiel d. 
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Camera  System 


9.  A  35  mm  SLR  remote  control  camera  equipped  with  a  flash  and 
self  winder  was  used  during  model  studies  of  Phases  I  and  II.  The 
camera  recorded  a  cross-sectional  view  of  the  test  embankment  shown  in 
the  mirror.  Additionally  all  tests  were  viewed  in  real  time  using  the 
closed  circuit  video  camera,  and  recorded  on  video  tape. 

Water  Depth  Measurement 


10.  In  the  rigid  model  tests,  the  depth  of  overtopping  was 
monitored  by  measuring  the  voltage  drop  across  two  copper  plates 
immersed  in  the  flow.  In  addition  the  water  depth  was  also  measured  in 
reference  to  the  grid  system  inscribed  on  the  plexiglass  sidewall,  both 
by  visual  observation  on  the  video  screen  during  the  test  and  by 
post-test  analysis  of  photographs  taken  in  flight.  The  two  methods 
were  in  close  agreement  in  all  tests.  In  the  erodable  model  tests, 
only  the  grid  system  was  used  for  depth  measurements. 
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Figure  1 1 - 1  Schematic  Of  10  g-ton  Geotechnical  Centrifuge 
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direction  of  rotation  in  the 
centrifuge 


Figure  1 1-3  Schematic  Of  Sample  Container  In  Its  Inflight 
Position,  Plan  View 
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figure  I 1-4  Schematic  Of  Sample  Container  In  Its  Inflight 
Position, Side  View 
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Figure  1 1 -5  (a)  Initial  Mounting  Position  Attempted 
(b)  Mounting  Position  Utilized 


PART  III:  TEST  MATERIAL  AND  PREPARATION 
OF  SAMPLE  EMBANKMENTS 

Phase  I  -  Rigid  Embankments 


11.  The  configuration  of  the  prototype  SLA  embankment  is  shown 
in  Figure  III-l.  This  configuration  was  simplified  so  that  the  crest 
was  flat  rather  than  having  the  compound  slope  shown  in  the  figure. 
Three  scale  models  of  this  embankment  were  prepared:  l/24th,  l/49th, 
and  1/7 1st  scale.  When  accelerated  at  24,  49,  and  71  g's,  respec¬ 
tively,  these  models  represent  the  prototype  SLA  embankment.  All 
models  were  constructed  of  aluminum  and  the  flow  surfaces  were  machined 
using  an  end  mill  to  create  a  smooth  surface.  As  indicated  previously 
due  to  space  limitations  only  a  portion  of  the  crest  and  the  downstream 
slope  were  modeled.  Additionally,  in  the  l/24th  model  no  toe  section 
was  modeled.  The  configurations  of  the  three  models  are  shown  in 
Figures  1 1 1- 2  through  1 1 1-4. 

12.  The  front  sides  of  the  embankments  were  painted  flat  black, 
while  the  backwall  of  the  sample  container  was  painted  pink,  to  provide 
contrast  with  the  overtopping  water  which  was  dyed  blue,  thus  facili¬ 
tating  measurement  of  the  flow  overtopping  depth  in  the  photographs. 

Phase  II  -  Erodable  Embankments 


Soil  Material 

13.  The  sandy  clay  soil  used  by  SLA  for  prototype  embankment 
construction  was  also  used  for  construction  of  the  erodable  embankment 
models.  Soil  properties  and  the  grain  size  curve  for  this  material  are 
shown  in  Table  1 1 1 - 1  and  Figure  1 1 1 - 5  respectively.  The  grain  size 
distribution  was  not  scaled  downward  for  two  reasons.  First,  the 
behavior  of  clay  soils  is  primarily  dependent  on  surface  forces,  rather 
than  gravity  induced  body  forces,  due  to  their  small  particle  size. 
Second,  scaling  of  the  grain  size  distribution  would  require  an 
elaborate  centrifugation  process  to  separate  the  soil  grains  for 
reconst ituti on. 
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Preparation  of  the  Test  Embankments 


14.  The  soil  obtained  from  the  field  was  oven  dried  at  110°C, 
and  crushed  in  a  mechanical  crusher  and  passed  through  a  No.  8  sieve  to 
eliminate  any  large  clods  or  pieces  of  gravel.  Sufficient  dry  material 
to  form  a  single  embankment  was  mixed  in  an  electric  mixer  with  the 
appropriate  quantity  of  water  to  obtain  an  optimum  water  content  of 
14%.  The  material  was  then  compacted  statically  using  a  hydraulic 
press  by  the  undercompaction  method  (Ladd,  1978)  in  0.5  inch  layers  in 
the  sample  container  to  a  dry  density  of  106  pcf  (90%  of  Modified 
Poctor  Maximum  Dry  Density)  at  a  moisture  content  of  14%  into  a 
rectangular  block  as  shown  in  Figure  1 1 1 -6 .  The  surface  of  each  layer 
was  thoroughly  scarified  prior  to  placing  the  subsequent  layers. 

Density  was  verified  by  weighing  the  sample  container  after  completion 
of  compaction.  Following  the  compaction,  the  plexiglass  front  wall  was 
removed  and  the  embankment  configuration  was  trimmed  using  a  sharp 
knife  and  template  to  the  same  configuration  used  in  the  three  rigid 
embankments  shown  in  Figures  1 1 1 -2  through  II 1-4.  In  addition  to  the 
models  with  the  same  geometry  as  the  rigid  embankments,  one  additional 
test  was  performed  to  model  the  entire  geometry  of  the  full  scale  SLA 
test  at  1/7 1st  scale.  The  single  1/7 1st  scale  test  embankment 
constructed  with  upstream  and  downstream  slopes  is  shown  in  Figure 

1 1 1-7. 

15.  The  sample  container  backwall  was  painted  pink  as  in  the 
case  of  the  rigid  embankments,  in  order  to  provide  contrast  with  the 
blue  overtopping  water  and  brown  soil. 

16.  The  compaction  procedure  used  for  embankment  preparation  was 
different  than  that  used  by  SLA  for  the  actual  embankment  sections. 
According  to  information  obtained  from  SLA  the  field  embankments  were 
compacted  directly  in  the  test  flume.  The  material  was  compacted  at  a 
moisture  content  of  18%  in  6  inch  thick  lifts  using  a  manual  jumping 
jack  type  compactor  to  95%  of  Standard  Proctor  Maximum  Dry  Density. 

Thus  while  model  and  prototype  densities  were  comparable,  compaction 
varied  in  lift  thickness,  moisture  content,  and  compaction  method. 
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TABLE  III-l 


PROPERTIES  OF  SANDY  CLAY  TEST  SOIL 

UNIFIED  SOIL  CLASSIFICATION  CL/ML  &  OL 

GRAIN  SIZE  DISTRIBUTION 
(100%  passing  #4) 

Percent  Sand  (+  #200)  37 

Percent  Si  It  50 

Percent  Clay  (-2pm)  13 

ATTERBERG  LIMITS 

Liquid  Limit  39 

Plastic  Limit  26 

Plasticity  Index  13 

COMPACTION  PROPERTIES 

(Modified  Proctor) 

Maximum  Dry  Density  118  pcf 

Optimum  Water  Content  14% 
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3.50 


Figure  1 1 1-2  1/24 th  Scale  Rigid  Embankment,  1.50"  In  Thickness 


18 


Scale  Rigid  Embankment,  0.74”  In  Thickness 


Figure  1 1 1 -5  SLA  Embankment  Soil  Grain  Size  Distribution 


21 


22 


st  Scale  Erodable  Full  Embankment  Section  Model 


PART  IV:  TEST  PROCEDURE 


Phase  I  -  Rigid  Embankment  Tests 


17.  Three  rigid  embankment  models  were  used  in  this  segment  of 
the  study  to  perform  a  "modeling  of  models"  study,  and  to  evaluate 
modeling  of  flow  quantities,  from  which  flow  velocities  can  be  readily 
calculated.  After  installation  of  an  embankment  model  in  the  test 
container,  this  assembly  was  placed  in  the  centrifuge  payload  basket  in 
the  configuration  shown  in  Figure  1 1 -4.  One  35  mm  camera  with  flash  in 
combination  with  the  video  television  camera  was  used  to  view  the  flow 
cross  section  through  the  mirror.  After  taking  a  picture  with  no 
overtopping  water,  the  centrifuge  was  spun  up  to  the  speed  necessary  to 
produce  the  desired  gravity  level  required  for  similarity  of  each 
model.  Once  at  that  speed  the  video  tape  recorder  was  turned  on  to 
record  the  experiments.  Then  the  water  pump  was  turned  on  to  introduce 
water  into  the  centrifuge.  The  centrifuge  speed  was  adjusted  as 
necessary  to  maintain  the  proper  g  level.  The  overtopping  depth  was 
determined  using  the  electrical  depth  meter  installed  in  the  back 
compartment  reservoir  and  the  grid  system  inscribed  on  the  plexiglass 
wall  of  the  test  container.  Photos  were  taken  by  remote  control,  and 
the  flow  quantity  was  measured  using  the  method  described  previously. 
The  overtopping  depth  was  then  increased  to  a  higher  value,  the 
centrifuge  speed  adjusted  as  necessary,  photos  taken,  and  the  flow 
quantity  measured.  After  sufficient  measurements  were  taken,  the  water 
pump  was  turned  off  and  the  centrifuge  was  slowly  brought  to  a  stop. 

Phase  II  -  Erodable  Embankment  Tests 


18.  Three  erodable  soil  embankment  models  were  used  in  this 
study  segment  to  evaluate  erosion  modeling,  again  using  the  modeling  of 
models  technique.  After  compaction  and  trimming  of  the  embankment 
model,  the  test  container  was  assembled  and  installed  in  the  centrifuge 
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using  the  same  configuration  utilized  in  Phase  I  tests.  Figure  1 1 -4. 
Both  the  single  35  mm  camera  and  video  camera  were  used  for 
experimental  recording.  After  taking  an  initial  photo  of  the 
embankment  prior  to  overtopping,  the  centrifuge  was  spun  up  to  the 
proper  speed  for  the  desired  gravity  level,  and  the  water  pump  was 
switched  on.  As  quickly  as  possible  the  flow  control  valve  was 
adjusted  to  achieve  overtopping  flow  depths  of  1.0  inch  in  the  l/24th, 
0.49  inches  in  the  l/49th,  and  0.34  inches  in  the  1/7 1st  scale  models. 
These  depths  all  model  a  2.0  ft.  overtopping  depth  in  the  SLA  prototype 
embankment. 

19.  Immediately  upon  reaching  the  desired  overtopping  depth, 
taken  as  zero  time,  a  stopwatch  was  started.  Pictures  were  taken 
during  overtopping  at  the  rate  of  about  one  per  minute  to  allow 
monitoring  and  measurement  of  erosion.  After  overtopping  was 
completed,  the  water  pump  was  turned  off,  the  centrifuge  was  slowly 
stopped,  and  the  overtopping  flow  measured  using  the  method  described 
previously. 
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PART  V:  EXPERIMENTAL  RESULTS 
Phase  I  -  Rigid  Embankment  Tests 


20.  In  this  phase  of  the  experimental  program  tests  on  rigid 
non-erodable  embankment  models  were  completed  to  investigate,  with  a 
modeling  of  models  analysis,  the  overtopping  depth  versus  flow  quantity 
relationship,  from  which  flow  velocity  can  easily  be  determined. 

21.  Overtopping  depths  were  measured  by  mounting  in-flight 
photograph  negatives  in  slide  mounts  and  projecting  them  at  actual 
scale  onto  a  sheet  of  paper.  The  embankment  and  overtopping  flow 
cross-sections  were  then  traced  on  the  paper,  and  measurements  then 
taken  directly  off  these  drawings.  Data  from  the  tests  used  in  the 
analysis  is  tabulated  in  Table  V-l,  and  corresponding  photographs  of 
the  overtopping  in  the  l/24th,  l/49th,  and  l/71st  scale  models  are 
presented  in  Figures  V-l,  V- 2,  and  V-3  respectively.  Overtopping 
depths  H  were  measured  in  the  reservoir  behind  the  embankment  models. 
Flow  quantities,  measured  by  the  method  previously  discussed,  were 
divided  in  half  as  it  was  assumed  that  one-half  the  flow  was  directed 
to  the  counterweight  basket.  This  adjusted  value  is  listed  as 
overtopping  flow  Q  in  Table  V-l.  The  widths  of  the  models  are 
presented  in  the  table.  As  shown  in  Figure  V-4,  overtopping  flow  in 
the  back  compartment  had  to  turn  90°  to  enter  the  embankment  flow 
section  and  it  was  observed  that  this  led  to  a  narrowing  of  the  flow 
width  similar  to  that  described  by  Chow  (1959)  for  contractions  in  a 
sharp  crested  weir.  Therefore,  the  channel  width  was  corrected  to  an 
effective  width  using  the  following  equation  proposed  by  Chow: 

L  =  W  -  0.1NH  (V-l) 

in  which  W  is  the  actual  embankment  width,  N  is  the  number  of 
contractions  (taken  as  1  in  this  study)  and  H  is  the  overtopping 
depth.  The  flow  intensity  q  (cfs/ft)  was  calculated  by: 
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q=£  (V-2) 

22.  In  order  to  evaluate  modeling  of  the  relationship  between 
flow  intensity  q  and  overtopping  depth  H,  the  data  were  plotted 
logarithmically  as  shown  in  Figure  V-5,  and  a  least-squares  analysis 
was  used  to  fit  curves  to  the  data.  The  resulting  curves  were  found  to 
fit  the  general  form: 

q  =  anH*  (V-3) 

where  6  =  1.5  and  an  i  s  a  coefficient  related  to  the  model  gravity 
level.  From  the  data  shown  in  Figure  V-5,  a2lt  =  16.0,  a49  =  35.0, 
d7i  =  44.7.  This  coefficient  an  is  governed  by  the  following  scaling 
relationship: 

ainiX  =  annX  (V-4) 

By  rearranging  terms,  and  recognizing  that  n^l  for  the  prototype 
situation, 

an  =  ain-X  (V-5) 

where  at  is  the  prototype  coefficient. 

n  =  (V- 6) 

9p 

where  gm  and  gp  are  the  model  and  prototype  accelerations 
respectively.  The  study  data  indicated  \  =  "0.97. 

23.  The  embankment  sections  tested  are  essentially  broad-crested 
weirs,  and  in  fluid  mechanics  flow  over  this  type  of  section  is  usually 
defined  as  (Chow,  1959): 

q  =  Cd  /Zg  H3/2  (V- 7) 

which  when  modified  for  varying  accelerations  becomes: 

Q  =  Cdn  /2grT  H3/2  (V-8) 
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where  C^p  is  an  experimentally  determined  coefficient  of  discharge 
dependent  upon  the  entrance  geometry  of  the  weir  and  on  n.  All  other 
terms  are  as  previously  defined.  Comparing  equations  V- 3  and  V- 8  and 
noting  that  the  value  of  $  found  from  the  experimental  data  was  indeed 
equal  to  that  proposed  by  Chow  (1959),  we  can  see  that 

a„  -  Cd„  /2gn  (V-9) 

Application  of  this  relationship  to  the  experimental  data  allowed  us  to 
calculate  C^ n  and  to  plot  it  as  shown  in  Figure  V-6,  from  which  the 
following  scaling  relationship  for  C^n  can  be  developed: 

cdl  n/  -  C dy  (V-10) 

and  rearranging  terms: 

Cdn  “  Cdl"'*  ‘V-U> 


where  Cdl  is  the  discharge  coefficient  value  of  the  prototype.  The 
model  data  indicated  y  =  -0.47  and  C^i  =  0.094.  If  equation  V - 1 1  is 
substituted  for  C^n  in  equation  V-8  the  relationship  for  flow  then 
becomes: 


q 


(V-12) 


and  substituting  experimentally  determined  values  for  C^i  and  y,  and 
for  g  the  flow  relationship  was  found  to  be: 

q  =  o.74  n°*97H3/2  (V-13) 

which  if  rounded  off  becomes, 

q  =  0.75  nH3/2  (V-14) 

This  equation  was  found  to  fit  the  experimental  data  shown  in  Figure 
V- 5  very  wel 1 . 
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Phase  II  -  Erodable  Embankment  Tests 


24.  In  the  second  phase  of  the  study  tests  were  completed  on 
models  of  the  SLA  prototype  embankments  to  develop  data  for  evaluation 
of  centrifuge  modeling  of  erosion  rates.  To  allow  for  development  of 
scaling  relationships  governing  erosion  rates,  the  modeling  of  models 
technique  was  used  in  this  phase  in  tests  on  l/24th,  l/49th,  and  1/7 Is t 
scale  models. 

25.  Inflight  photographs  of  each  test  were  taken,  and  one 
representative  sequence  from  each  size  of  model  is  presented  in  Figures 
V-7  to  V- 9.  The  photos  show  that  in  each  test  erosion  started  at  the 
toe  of  the  embankment  slope  and  progressed  upward  toward  the  crest, 
creating  a  near-vertical  erosion  scarp.  Minor  surface  erosion  was  also 
observed  to  occur  away  from  the  toe  area,  which  is  thought  to  be  due  to 
poor  compaction  layer  contact  and  minor  surface  loosening  of  material 
during  trimming  of  the  embankment. 

26.  In  order  to  measure  the  eroded  soil  volumes,  the  negative  of 
inflight  photos  were  projected,  at  full  scale,  on  to  sheets  of  paper, 
and  the  eroded  shape  was  then  outlined.  The  original  non-eroded  shape 
was  also  traced  onto  each  sheet  to  facilitate  measurements.  Because 
the  l/24th  scale  embankment  did  not  include  a  toe,  its  area  at  appro¬ 
priate  scale  was  traced  on  to  drawings  from  the  l/49th  and  1/7 1st  scale 
embankments  as  shown  in  Figure  V-10,  to  define  the  "erodable  volume" 
for  all  tests  in  order  to  facilitate  comparison  of  data.  The  eroded 
areas  A  were  measured  with  a  polar  planimeter  and  erosion  volume  V  was 
assumed  to  be  equal  to  the  product  of  erosion  area  A  and  the  embankment 
thickness.  The  eroded  volumes  at  any  time  of  erosion  were  then  normal¬ 
ized  by  dividing  by  the  total  "erodable"  volume,  to  calculate  percent 
of  embankment  eroded.  Calculated  data  is  tabulated  in  Table  V-2  and 
plotted  in  Figure  V- 1 1 .  In  the  l/49st  scale  tests,  data  from  all  three 
tests  exhibit  good  agreement,  while  in  the  1 / 2 4t h  and  1/7 1th  scale 
models  two  tests  agree  well  while  the  third  test  deviated.  It  is  not 
known  why  this  occurred;  however,  variations  in  compaction  from  test  to 
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test  are  suspected.  Data  from  the  two  tests  that  deviated,  24-2  and 
71-2,  were  not  included  in  development  of  modeling  relationships  for 
erosion.  The  erosion  curves  shown  in  Figure  V-ll  show  behavior  in 
which  erosion  was  progressing  at  a  fairly  steady  rate,  then  suddenly 
increased  due  to  a  large  loss  of  material  and  then  resumed  at 
essentially  the  original  or  similar  rate.  For  example,  in  test  24-1, 
erosion  was  progressing  up  to  16  minutes  at  a  steady  rate  and  in  the 
next  2  minutes  a  large  volume  of  soil  eroded.  After  19  minutes  erosion 
rate  resumed  at  level  similar  to  that  before  16  minutes.  Furthermore, 
in  all  tests  there  seemed  to  be  a  primary  erosion  rate  up  to  4  minutes, 
then  followed  by  an  essentially  linear  segment,  if  any  sudden  loss  of 
soil  is  neglected.  The  slope  of  these  linear  portions  are  defined  as 
the  erosion  rates  and  were  used  for  the  modeling  of  models  analysis. 

27.  The  erosion  depths  at  six  points  on  the  embankment  were  also 
measured  at  the  locations  indicated  in  Figure  V-12,  and  were  converted 
to  prototype  scale  for  analysis.  The  data  from  these  measurements  is 
presented  in  Figures  V-13  through  V-18.  The  variability  of  the  data 
appears  to  be  fairly  random,  which  tends  to  indicate  variations  in 
compaction. 

28.  This  variability  also  indicates  that  slope  erosion  is  a 
process  with  a  certain  degree  of  randomness.  The  data  shows  that 
during  the  erosion  process  material  is  removed  by  both  fairly  steady 
surface  erosion,  and  by  discrete  events  in  which  large  chunks  of 
material  are  removed  quite  rapidly.  The  data  contained  in  Table  V-3 
represent  erosion  rates  for  the  model  embankments. 

29.  The  average  rate  for  each  model  scale  was  calculated  as 
shown  in  Table  V-3  and  the  following  relationship  between  model  erosion 
rate  Em  i n  %  vol/min  and  scaling  factor  was  obtained. 

E  =  0.0077  n1*38  (V-15) 

m 

where  model  erosion  rate  is  defined  as: 

E  =  V^°^d/Vori9inal  .  V^o  .  %  vol/m1n  (y-16) 

m  t  t 

The  prototype  erosion  rate  (n=l)  was  determined  to  be: 
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Ep  =  0.0077  %  vol/min 


(V- 17 ) 


and  thus  the  scaling  relation  for  erosion  rate  was  found  to  be: 

c  n  nrt->7  1.38 

E  0.0077  n  .  .q 

P  P  “ 1  •  jO 

r  —  r  =  n 


E  0.0077  n 
m 


1.38 


Substituting  equation  (V-15)  into  equation  (V-18): 


m 


(V  /V  ) 
v  e  o'p 

(V  /V  ) 
v  e'  o'm 


m  -1.38 
—  =  n 


(Vp  (Vm  1.38 
•  n 


m 


<Vm  (Vo)p 


(V-18) 


( V— 19 ) 


the  scaling  relationship  for  time  in  the  erosion  process  can  be 
determined  from  the  experimental  data  as: 
t 


(V-20) 


The  scaling  relationship  for  volume  is  1  to  n3  so  equation  (V-20)  can 
be  transformed  to: 


1.38 


=  n 


1.38 


(V-21) 


m 


Full  Embankment  Section  Test 


30.  One  1/7 1st  model  of  the  full  embankment  section  being  tested 
by  SLA  was  modeled  in  the  centrifuge  to  allow  direct  comparison  of 
erosion  data  with  the  prototype  data.  The  prototype  overtopping  depth 
used  was  2.0  ft.,  along  with  a  zero  tailwater  condition.  The  erosion 
data  provided  by  SLA  for  comparison  was  from  a  test  in  which  the  total 
duration  of  overtopping  was  17.5  hours.  Inflight  photographs  of  the 
model  are  shown  in  Figure  V- 19,  and  downstream  slope  cross-sections 
from  the  prototype  and  model  for  a  range  of  overtopping  times  are 
presented  in  Figure  V-20.  Erosion  of  the  crest  was  minimal  in  both  the 
model  and  prototype.  Cross-sections  were  reviewed  to  determine  those 
for  which  model  prototype  agreement  was  best.  As  shown  in  Figure  V-21 
prototype  erosion  at  2.5  hours  (150  min.)  agreed  very  well  with  model 
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cross-sections  at  1  and  3  minutes.  Evaluation  of  the  data  to  determine 
the  scaling  relationship  for  time  yields  the  following  relationships: 

—  =  n*'*8  (Model  time  =  1  minute)  (V-22) 

—  =  n8'^  (Model  time  =  3  minutes)  (V-23) 

tm 

Also  shown  in  Figure  V-21  are  erosion  cross-sections  for  model  time  of 
8  minutes  and  prototype  time  of  13.5  hours.  As  indicated  by  the  data, 
agreement  is  not  as  good  as  at  the  earlier  times,  although  the  shapes 
of  the  eroded  downstream  slopes  are  very  similar.  This  later  data 
indicates  the  following  scaling  relationships: 


(V-24) 


31.  The  relationships  shown  in  equations  V-22  through  V-24  do  not 
agree  with  that  determined  in  the  modeling  of  models  study  indicated  by 
equation  V-21.  If  the  scaling  relation  of  equation  V-21  is  used,  then 
the  model  time  corresponding  to  13.5  hrs  (810  minutes)  of  prototype 
time  is  calculated  to  be  2.26  minutes.  No  data  was  available  at  2.26 
minutes;  however,  a  profile  at  2.0  minutes  is  shown  in  Figure  V-21  for 
comparison.  It  can  be  seen  that  considerable  differences  are  found. 
However,  if  variations  in  density  and  soil  structure  and  thus  soil 
erodability  between  the  model  and  prototype  occurred,  this  would  result 
in  errors  in  the  determination  of  time  scaling  relationships  of 
equations  V-22  through  V-24.  As  discussed  previously  while  model 
densities  were  comparible  to  those  attained  in  the  prototype,  the 
compaction  procedures  varied  in  lift  thickness,  compaction  moisture 
content,  and  compaction  method.  These  differences  would  likely  result 
in  variations  in  soil  structure  between  the  model  and  prototype 
embankments.  However,  the  degree  of  this  variation  cannot  be 
determined  with  the  data  available. 
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TABLE  V-l 


RIGID  EMBANKMENT  TESTS  OVERTOPPING  DATA 


Model 

Scale 

Factor 

n 

Channel 

Width 

W 

(ft) 

Corrected 

Channel 

Width 

L 

(ft) 

Overtoppi ng 
Depth 

H 

(ft) 

Overtoppi ng 

FI  ow 

Q 

(cfs) 

Unit 

Overtopping 

Flow 

q 

(cfs/ft) 

24 

.125 

.118 

.067 

.031 

.26 

.119 

.058 

.026 

.22 

.119 

.060 

.027 

.22 

.117 

.079 

.041 

.35 

.116 

.090 

.049 

.42 

.115 

.098 

.055 

.47 

49 

.062 

.058 

.036 

.014 

.25 

.057 

.044 

.019 

.32 

.056 

.052 

.024 

.43 

.055 

.061 

.029 

.52 

.054 

.074 

.038 

.70 

.053 

.086 

.049 

.92 

.052 

.098 

.054 

1.1 

71 

.042 

.039 

.033 

.011 

.28 

.037 

.049 

.018 

.48 

.035 

.069 

.027 

.76 

.034 

.078 

.034 

1.0 

.032 

.098 

.038 

1.2 
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TABLE  V-2 


EROOABLE  EMBANKMENT  DATA 


Test 

1 

Test 

2 

Test 

3 

Model 

Overtoppi ng 

Eroded 

Overtopping 

Eroded 

Overtopping 

Eroded 

Scale 

Time 

Vol  ume 

Time 

Volume 

Ti  me 

Vol  ume 

Factor 

t 

V 

t 

V 

t 

V 

n 

(min) 

(%) 

(mi  n) 

(%) 

(mi  n) 

(%) 

24 

1.25 

3.7 

1.00 

2.2 

1.50 

0.9 

2.33 

6.8 

2.00 

5.5 

2.17 

1.6 

4.58 

14.6 

4.00 

21.6 

4.17 

7.8 

8.75 

23.0 

8.00 

54.6 

8.00 

13.7 

16.50 

28.1 

15.30 

81.5 

16.00 

20.8 

19.75 

48.1 

23.00 

88.7 

26.00 

26.0 

25.00 

53.9 

36.00 

28.8 

30.00 

55.4 

46.00 

35.2 

36.00 

59.1 

49 

1.00 

4.5 

0.83 

1.5 

0.83 

0.9 

1.92 

6.5 

2.10 

4.4 

2.17 

3.7 

4.17 

7.9 

4.08 

5.6 

4.17 

8.2 

8.17 

9.9 

8.00 

9.3 

8.00 

33.2 

17.17 

12.4 

16.83 

11.0 

11.75 

38.7 

24.17 

47.4 

24.33 

19.6 

17.33 

40.3 

30.17 

47.8 

30.50 

25.8 

71 

1.00 

4.9 

1.00 

0.8 

1.00 

9.9 

2.00 

4.9 

2.00 

4.8 

2.00 

11.9 

4.00 

6.5 

4.00 

5.9 

4.00 

15.9 

8.50 

28.3 

8.50 

26.2 

8.00 

32.9 

16.00 

44.0 

11.00 

48.1 

13.00 

47.1 

19.00 

49.2 

13.00 

68.1 

16.00 

49.5 

15.50 

74.5 
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TABLE  V-3 


EROSION  RATE  DATA 


Model 

Scale 

Ti  me 

Range 

Erosion 

Rate 

Average 

Erosion 

Rate 

Factor 

n 

AT 

(mi  n) 

• 

E 

(  itvol/min) 

^ave 

(%vol/mi n) 

24 

19.75-36.00 

0.64 

0.62 

4.17-36.00 

0.60 

49 

4.17-30.17 

1.74 

1.57 

4.08-30.50 

0.73 

4.17-17.33 

2.24 

71 

4.00-19.00 

2.73 

2.79 

4.00-16.00 

2.85 
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Figure  V-l  l/24th  Scale  Rigid  Embankment,  (a)  H=0',  (b)  H=0 . 058 1 ,  (c)  H=0.060',  (d)  H=0.067 
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9th  Scale  Rigid  Embankment,  (a)  H=0',  (b)  H=0 . 036 1 ,  (c)  H=0 . 044 ' ,  (d)  H=0 .052 


Scale  Rigid  Embankment,  (e)  H=0.061’,  (f)  H=0.074',  (g)  H=0.086',  (h)  H=0.098 
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Figure  V-3  1/ 7 Is t  Scale  Rigid  Embankment,  (a)  H=0 1 ,  (b)  H=0 . 033 ' ,  (c)  H-0.049',  (d)  H=0.069 


Figure  V-4  Contraction  Of  Flow  At  Entrance  To  Embankment 
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Figure  V-5  Un 


n 


Figure  V-6  Discharge  Coefficient-  Versus  Model  Scale  Factor-  n 
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Figure  V- 7  1 / 24th  Scale  Lrodable  Embankment,  (e)  t=8.75  min,  (f)  t = 1 6 . 5  min,  (q)  t=  19.75,  (h)  t  =  25  min 


45 


w 


46 


Figure  V-8  l/49th  Scale  Erodable  Embankment,  (a)  t=0  min,  (b)  t- 1 .00  min,  (c)  t-1.92  min, 
(d)  t=4 .17  min 
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Figure  V-8  l/49th  Scale  Erodable  Embankment,  (e)  t=8,17  min,  (f)  t= 17 . 17  min,  (g)  t=24 . 17  min, 
(h)  t=30.17  min 
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Figure  V-9  1/7 1 s t  Scale  Erodable  Embankment,  (a)  t=0  min,  (b)  t=1.00  min,  (c)  t=2.00  min, 
(d)  t=4.00  min 
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Figure  V- 9  l/71st  Scale  Erodable  Embankment,  (e)  t=8.50  min,  (f)  t=16.00  min,  (g)  t=19.00  min 


Erodable  Volume 
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Figure  V- 10  Erodable  Volume  From  l/24th  Scale  Embankment  Superimposed  On  l/49th 
And  1/7 Is t  Scale  Embankment  Cross-Sections. 
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re  V- 13  Erosion  Depth  Versus  Time 


Erosion  Depth  (in) 


Erosion 


t  Point  3 


Depth  Versus  Time 
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100.0  r 
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Figure  V- 17  Erosion  Depth  Versus  Time  At  Point  5 
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Figure  V-19  1/7 Is t  Scale  Erodable  Full  Embankment 
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Figure  V- 19  l/71st  Scale  Erodable  Full  Embankment  Section:  (e)  t=9.0  min. 


SLA  Prototype  at  2.5,  5.5,  8.5,13.5,  15.5,  and  17.5  hours  of 
Overtopping  Flow 


1/7J  rCale  Model  at  ll0’  3-°*  7-2>  s-0’  9-°*  10.0,  14.5,  and 
19.5  hours  of  Overtopping  Flow  (Prototype  time) 

Figure  V-20  Prototype  and  Full  Embankment  Model  Downstream  Slope 

Erosion  Profiles  (Scale:  1.0  inch  =  6.0  feet,  prototype) 


-61- 


1/71 st  Scale  Model  at  1.0  and  3.0  minutes  of  overtopping 
SLA  Prototype  at  2.5  hours  of  overtopping 


Figure  V-21  Prototype  and  Full  Embankment  Model  Downstream  Slope 

Erosion  Profiles  (Scale:  1.0  inch  =  6.0  feet,  prototype) 
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PART  VI:  CONCLUSIONS 


Phase  I  -  Rigid  Embankment  Tests 


32.  In  order  to  show  the  applicability  of  the  centrifuge  to 
modeling  of  hydraulic  parameters  related  to  open-channel  flow  a 
modeling  of  models  analysis  was  completed  on  a  series  of  rigid 
embankments.  The  results  of  the  analysis  indicated  that  flow 
quantities  q  modeled  quite  well  in  the  centrifuge  and  the  following 
relationship  was  determined: 

q  =  0.75  n  h3/2  (V-14) 

This  equation  agrees  with  that  presented  by  Chow  (1959)  for 
broad-creased  weirs,  with  the  exception  that  experimental  discharge 
coefficient,  0.75  falls  below  the  range  of  2.46  to  3.47  presented  by 
Chow.  It  is  felt  that  this  difference  is  due  to  entrance  conditions 
used  in  the  centrifuge  models  necessitated  by  the  limited  space 
available.  These  conditions  brought  about  a  different  ratio  of  depth 
of  upstream  pool  to  embankment  height,  and  required  that  overtopping 
flow  had  to  turn  a  corner  just  prior  to  entering  the  flow  channel. 

Both  of  these  altered  states  would  change  the  flow  regime  upon  which 
Chow's  values  are  based. 

Phase  II  -  Erodable  Embankment  Tests 


33.  The  SLA  prototype  embankment  section  was  simulated  in  a 
modeling  of  models  scheme.  The  same  sandy  clay  soil  was  used  to 
construct  models  at  l/24th,  l/49th,  and  l/71st  scale  which  were  then 
tested  in  triplicates,  for  a  total  of  nine  tests.  In  all  of  these 
models  erosion  started  at  the  toe  and  progressed  upwards  towards  the 
crest,  a  mechanism  like  that  observed  in  the  feasibility  study  tests 
(Ko  et  al.,  1984). 

34.  Comparison  of  erosion  depth  data  at  various  points  along  the 
model  slopes,  indicated  a  high  degree  of  variability  from  model  to 


model,  which  was  attributed  to  minor  variations  in  density  in  the 
models  due  to  preparation  techniques  and  to  the  random  nature  of  the 
erosion  process.  When  data  on  percentage  of  embankment  volume  eroded 
was  compared,  variability  of  data  was  reduced. 

35.  The  modeling  of  models  analysis  completed  indicated  that 
erosion  rates  modeled  well,  and  the  following  relationship  for  model 
erosion  rate  in  %  vol/min.  was  determined: 

E  =  0.0077  n1*38  (V-15) 

m 

The  time  scaling  relationship  for  erosion  was  thus  found  to  be: 


m 


36.  One  1/7 1st  scale  full  embankment  section  model  was  tested 
and  the  erosion  data  was  compared  to  prototype  data  provided  by  SLA  to 
evaluate  the  time  scaling  relationship.  This  comparison  yielded 
scaling  factors  of  time  of  n8*92  to  n^*^8,  which  did  not  agree  with 
equation  V-21.  This  was  likely  due  in  part  to  variations  in  the 
structure  of  the  soil  in  the  model  from  that  in  the  prototype  as  a 
result  of  variations  in  compaction  technique. 
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